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In the derivation of low-energy effective models for solids targeting the bands near the Fermi level, 
the constrained random phase approximation (cRPA) has become an appreciated tool to compute 
the effective interactions. The Wick-ordered constrained functional renormalization group (cfRG) 
generalizes the cRPA approach by including all interaction channels in an unbiased way. Here 
we present applications of the cfRG to two simple multi-band systems and compare the resulting 
effective interactions to the cRPA. First we consider a multiband model for monolayer graphene, 
where we integrate out the cr-bands to get an effective theory for 7r-bands. It turns out that terms 
beyond cRPA are strongly suppressed by the different 2 ;j/-plane reflection symmetry of the bands. In 
our model the cfRG-corrections to cRPA become visible when one disturbs this symmetry difference 
slightly, however without qualitative changes. This study shows that the embedding or layering 
of two-dimensional electronic systems can alter the effective interaction parameters beyond what 
is expected from screening considerations. The second example is a one-dimensional model for a 
diatomic system reminiscent of a GuO chain, where we consider an effective theory for Cu 3d-like 
orbitals. Here the fRG data shows relevant and qualitative corrections compared to the cRPA 
results. We argue that the new interaction terms affect the magnetic properties of the low-energy 
model. 


I. INTRODUCTION 


The study of correlation and ordering effects in solids on 
lower energy scales is very often performed by using low- 
energy effective models (LEEM). In formulating such a 
model, much of the material-specific information has to 
be cast into the parameters of the LEEM. A standard ex¬ 
ample for this approach is the one-band Hubbard model. 
In its derivation (see e.g. Ref. [l|), a single band at the 
Eermi level, held responsible for correlation effects like 
magnetism and unconventional superconductivity, is cho¬ 
sen as single-particle basis for the LEEM. The hopping 
parameters of the Hubbard model encode the material- 
specific extent and energetics of Wannier states arising 
from the electronic and crystal structure. Similarly, the 
interaction parameter of the Hubbard model, U, depends 
on the electronic structure via Coulomb integrals of the 
Wannier basis, and additional screening by other pro¬ 
cesses, as will be detailed further below. Once the Hub¬ 
bard model is formulated, many-body methods capable 
of capturing electronic correlations to sufficient detail can 
be used to solve for the low-energy properties. Usually, 
this last step is already a complicated endeavor, and for 
a long time, most researchers have focussed on getting 
the qualitative features of the model physics righL e.g. 
to understand the Mott transition (see, e.g. Ref. In 
recent time however, quantitative aspects have attracted 
more interest again. The combination of ab-initio meth¬ 
ods with dynamical mean-field theory has become a lively 
field of activities (see, e.g. Ref. i and 0) ■ In the field of 
high-temperature superconductivity, various groups have 
tried to understand the systematics of the T^-differences 
between various high-Tc cuprates^^— or iron arsenidesi^d^ 


Likewise, in the field of graphene, an understanding of 
the ground state of a single graphene sheet as well as of 
its layered cousins requires a more precise understanding 
what the actual model parameters areJA— 

Regarding the effective model parameters, there may be 
a different degree of controlled knowledge about single¬ 
particle parameters, like hopping amplitudes, compared 
to what is usually known about the interaction param¬ 
eters. The single-particle parameters can, at least in a 
dressed form, be estimated to large extent from exper¬ 
iments like angle-resolved photoemission. De-Haas-van- 
Alphen and angular-dependent magneto-resistance oscil¬ 
lation techniques are additional established complemen¬ 
tary tools to study the renormalized single-particle prop¬ 
erties. In contrast with this, for the interaction param¬ 
eters, there are at best indirect probes to their renor¬ 
malized values, and even then the correspondence to pa¬ 
rameters of a chosen LEEM may involve additional the¬ 
oretical steps. Therefore, this part of the LEEM de¬ 
scription is much less certain, may however be of cru¬ 
cial importance to understand material-specific inter¬ 
action effects in quantitative form. Fortunately, there 
have been at least two useful theoretical approaches 
that try to compute the relevant interaction parameters 
starting from the ab-initio electronic structure, beyond 
the direct calculation of Coulomb integrals in a local¬ 
ized basis set. These are the constrained local density 
approximation ; and the recently strongly used con¬ 
strained random phase approximation (cRPA^i^^— ). In 
the latter approach, the Coulomb integrals in the lo¬ 
calized basis are altered by dielectric functions arising 
the additional screening due to the bands not included 
in the LEEM. This screening is obtained from comput¬ 
ing the particle-hole polarizability of the band structure 
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constrained to those processes that do not happen ex¬ 
clusively within the target bands of the LEEM, as those 
will be dealt with in the solution of the LEEM and should 
not be treated twice. As is clear from the name cRPA, 
this approach is an approximation with caveats that are 
shared by many other RPA approaches. It is clear how¬ 
ever that the cRPA captures important physical ingredi¬ 
ents and it has been a very welcomed innovation over the 
last years. 

Here we ask the question if one can improve on the cRPA 
in a systematic way, in order to get even more quantita¬ 
tive precision. From other instances in physics where 
RPA and related perturbation expansions are used, it 
is known that RPA in the charge channel with a po¬ 
larizability as result can overlook important parts of 
physics. Moreover, it is not clear to us which small 
parameter warrants that an RPA summation for the 
LEEM interactions is necessarily close to the full an¬ 
swer. In fact, that perturbation channels different from 
the particle-hole RPA terms lead to important renor¬ 
malizations has been well known in condensed matter 
physics for a long time. Anderson and More k^i^° showed 
that the particle-particle channel causes a strong suppres¬ 
sion of the electron-electron repulsion at energy scales 
above the phonon frequencies, opening the chance that 
weak phonon-induced interactions lead to Cooper pair¬ 
ing. Around the same time, Kanamori^i employed for¬ 
mally related calculations again exploiting the particle- 
particle channel in order to compute effective interac¬ 
tions for magnetic ordering. Hence, it is not surpris¬ 
ing that the particle-particle channel will also affect the 
low-energy effective interactions in a multi band system. 
A quantitative difference may be that in the multiband 
case treated in this paper, the ’lower cutoff’ at which the 
renormalizations of interest are computed is still of the 
order of electron volts, and not much below as, e.g., in 
the Cooper pairing problem. Hence the logarithmic con¬ 
tributions in the particle-particle channel do not come to 
full strength, and the corrections are not as strong. Yet, 
as we will show below, they can be significant. Moreover 
similar things with noticeable magnitude happen in the 
’magnetic’ crossed particle-hole as well. 

In an earlier publication, one of us laid out a constrained 
functional renormalization group (cfRG) formalism that 
goes beyond the cRPA by including the other channels 
mentioned above plus vertex corrections, but can be re¬ 
duced to cRPA by dropping all but one one-loop dia¬ 
gram. In this earlier work, only very simplistic models 
were considered to show that the more general scheme 
can indeed lead to different answers. Here we use the 
cfRG in two more realistic model settings in order to 
show under which circumstances differences in the LEEM 
interactions can occur. The two models are 1) a simple 
8-band model for a graphene sheet, that in addition to 
the two TT-bonding orbitals per unit cell that give rise to 
Dirac spectrum also contains the six sp^-cr-bonding or¬ 
bitals per unit cell; and 2) an one-dimensional dps-model 
that mimics copper oxides or similar materials. Here the 


hopping between d-orbitals on Cu sites goes through p- 
orbitals on O sites. The main findings of this work are 
that in model 1) the cRPA gives a rather faithful picture 
of the LEEM interaction. Here the different symmetry 
of targets (tt-) bands and the (ct-) bands that are inte¬ 
grated out is the key factor for the absence of major cor¬ 
rections. Visible albeit limited corrections can however 
be provoked by breaking the symmetry of the graphene 
plane with respect to inversion of the perpendicular co¬ 
ordinate, i.e. by a substrate or a perpendicular electric 
field. In model 2) no such symmetry protection of the 
cRPA exists and we find rather marked deviations in the 
cFRG. Notably, the pairing and magnetic particle-hole 
channel cause distinct new features in the LEEM inter¬ 
actions, and the down-screening of the density-density 
terms is actually weakened compared to the cRPA result. 
We also try to understand what physical consequences 
arise from the extra terms in the effective interaction. 
Here, the main finding is an enhanced tendency of the 
LEEM to develop a large-spin ground state. This shows 
that the LEEM beyond the cRPA might even lead to new 
qualitative physical effects. 

This paper is organized in the following way: In SecHIlwe 
give a brief introduction to the constrained fRG method 
and discuss how the frequency structure of the vertex 
function can be interpretated in terms of density-density 
and spin-spin interactions. In Sec. IHIl we describe our re¬ 
sults for the graphene and the Cu-O-chain model. SecITVl 
contains the conclusions and a brief outlook. 


II. FORMALISM 

Here we review the Wick-ordered fRG formalism that, as 
argued in Ref. [12, can be used to compete the LEEM 
interactions beyond cRPA. 

A. Wick-ordered fRG 

The basic setup of what follows is an electronic struc¬ 
ture that can be divided in two sets of bands. The goal 
is to integrate out one set of bands (called ’high-energy 
bands’, although an energetic hierarchy is not essential 
for the formalism) in order to compute an effective de¬ 
scription for the remaining bands that are usually close 
to the Fermi level and that shall be called ’target bands’. 
In order to integrate out the high-energy bands we use 
the Wick-ordered fRG approach, described in Ref. [12, 
which lives in the functional integral formalism that de¬ 
scribes the band electrons as anti-commuting Grassmann 
fields. The quantum numbers carried by these fields are 
then a fermionic Matsubara frequency , a wave vector 
k in the first Brillouin zone of the lattice, the band-index 
b and a spin projection s € {f, i}. The band energies are 
then given by eb{k). We assume spin rotational invari¬ 
ance, i.e. the band are spin-degenerate. To parameterize 
the two-particle vertex function (ie. ’the interaction’) we 
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use bosonic ’Mandelstam’ Matsubara frequencies, which 
are related to the fermionic frequencies of the incoming 
(1 and 2) and outgoing (1’ and 2’) particles by the total 
incoming frequency ivi = za;i+iw2, and the two transfers 
ii>2 = iiov — iuii and = iui2 — ioJv. The SU( 2 )-spin- 
rotational invariance is reflected in our notation for the 
vertex function by the requirement that particle 1 and 1 ’ 
have the same spin projection si = ±i, and 2 and 2’ the 
same S2 = ±5- 

The fRG flow has the goal of integrating out the high- 
energy bands. It is organized in terms of Wick-ordered 
correlation function a^^i^^ and employs a generalized cut¬ 
off A that, upon evolving from 1 to 0 , removes the modes 
of the high-energy band from the theory and incorporates 
their effect in the remaining low-energy theoryi^ Com¬ 
pared to the usual amputated correlation functions that 
one obtains in the effective action by integrating out high- 
energy bands straightforwardly, the Wick-ordered two- 
particle correlation function is used as interaction of the 
LEEM. The Wick-ordered correlation function of n-th or¬ 
der contains in addition contributions from tree diagrams 
with n-|-2m external legs built with high-energy propaga¬ 


tors and one-particle-irreducible (IPI) vertex functions, 
where 2m legs have been contracted with m low-energy 
propagators. This resummation is crucial for incorpo¬ 
rating the ’mixed’ loop diagrams with one line in the 
high-energy and the second line in the low-energy sector. 
For more concrete formulae regarding the connection be¬ 
tween the different vertex functions, see the paper or 
textbooks like [H and [2^ 

For the integration over the modes of the high-energy 
bands the flow equations for the Wick-ordered ver¬ 
tices can be written down. As we focus in this 
work on the effective interactions, we will neglect 
self-energy corrections, although they may as well be 
of importance. Then fRG flow equations for the 
Wick-ordered two-particle interaction or coupling func¬ 
tion V^{ivi,iv2,ii'3',ki>,k2',ki-,bi>,b2’,bi,b2) are given 
by 23 ^ 

+ ‘j’crph (1) 

with 


^pp{ivi,iv2,iv3;ki>,k2',ki-,bi',b2',bi,b2) =-r'^ ^ — D^{iuj,k,b)D^(ivi - iuj^ky + k2' - k,b') 


b,b',k 


X _— fw, _^ — iio\ fci', ^2', k; by, & 2 ', b, b') 

X V^{ivi,iu} — iv^ _, _^ — iu)-, k, ky k 2 > — k,ky, b, b', bi, 62) 


^dph(**^i>**^ 2 ,W 3 ;*i',A: 2 ',A:i;&i/, 62 ',&i,& 2 ) ^ ^ D^{iuj,k,b)D^(iv2 + it^jky - ki + k,b') 


( 2 ) 


b.b',k 


_+ ioj — iv^ _; ky ^ fc, ky 6, 6i, h') 

xV^{iuj + iv+++, iv2, — fw; ky — ki + k, y, k; b', b2',b, 62) 

-I- V^{iu +^— -I- ioj, iv 2 , ioJ — -; ky, k, ky, by,b, bi, b') 

X V^{ioj + iv+++,iuj — iv+-+, —iv 2 \ ky —ki+k, fe', fe; ^2', &2, b) 

+ + iv+3 ,iui — iv^ _ k, ky, ky b, by,bi, b') 

X + iv+++,iv 2 ,'iv+-+ — ioj; ky — ki +k, k2',k-, b', b2',b, 62)] 


( 3 ) 


1 d [ ___ 

^ap\i{^'^i^^'^2,W3',ky,k2yky,by,b2',bi,b2)=-'^ ^ ^ D^{iuj,k,b)D^(iv3 + iuj,k + ky - ki,b') 


b,b',k 
rlk! 


X + iiy^ _ I-, iuj — iv^ _, W3; k, ky, ky b, by, bi, b') 

X V^iiuj + W+++, _— iuj, 11x3-, ky,k + y — ky by, b', b, 62) ( 4 ) 


Here we have defined _^ = \(ivi — iv2 +1^3) (w+++, 

_, w_|_are defined in an analogue way). Restrict¬ 
ing the flow to the underlined part of the direct particle- 
hole channel is equivalent to the cRPA solution.— In 


Fig. [T] we show a diagrammatic representation of the 
flow equations 
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The propagator is given by 


D^{b) 


(1 — A)C/o if & £ high energy bands 
Go if & £ low energy bands 


( 5 ) 




FIG. 1. (color online) Diagrammatic representation of the 
£RG flow equations, in the truncation that drops the six-point 
vertex. The dashed line carries a differentiated cutoff function 
and is only nonzero in the high-energy bands. The second loop 
line is a free propagator that can be in both high-energy and 
low-energy bands. Restricting the flow to the first diagram of 
the direct particle-hole channel in the second row reproduces 
the cRPA approximation. 


J 


As mentioned above, the flow equations are integrated 
from A = 0 to A = 1, in order that in the end of the flow 
the high energy bands are integrated out. Here two dif¬ 
ferences to usual momentum-shell RG schemes should be 
noted: First, we use a ’flat’ cutoff that simply switches off 
the high-energy bands in the free part of action irrespec¬ 
tive of their precise energy, and second, as characteristic 
feature of the Wick-flow, the second line in the one-loop 
terms is a low-energy propagator complementary to what 
has already been integrated out. Using a flat cutoff sim¬ 
plifies the numerical treatment significantly as long as no 
self-energy corrections are included. Then the diagrams 
on the right hand have to be commuted only once, and 
just get multiplied by —(1 — A) during the flow. 

It is actually possible to treat an interaction vertex that 
depends on three frequencies (see, e.g., Ref. but this 
leads to a larger numerical effort. In this study we choose 
to simplify the flow equations in the frequency domain 
and apply a Karrasch-Husemann decompositionj^i^ 
Here, the change of the interaction with A is written 
as a sum of three functions that each depend on only 
one of the the three bosonic frequencies mentioned be¬ 
fore, while the other two frequencies are kept at zero in 
these functions. Then the full vertex that still depends 
on three frequencies but can be written as a sum of three 
terms that each depend on one frequency. We conclude 
that the decomposition certainly suffices to study effects 
beyond RPA, but that at stronger couplings, even more 
frequency structures might emergei^i^ In any case, this 
approximation leads to the simplified flow equations 


--- I d T 

^pp{ivi;ki',k 2 ',ki;bv,b 2 ',bi,b 2 ) ^ ^ D^{iuj,k,b)D^(ivi - iuj,kv + k 2 '- k,b') 

b,b',k 

X V^{ivi,ivi /2 — iuj, iv\l 2 — ioj; fciq ^ 2 ', k; by, & 2 ', b, b') 

X V^{ivi,ioj — Wi/ 2 , wi /2 — fw; k, ky + ^ 2 ' — k,ki, b, h', bi, 62) 

^d-p\ii^'^ 2 ]ky,k 2 ',ki,by,b 2 ',byb 2 )=^Yl X! ^ D^{iuj,k,b)D^{iv 2 +ii^,ky - ki+k,b') 

b,b',k 

[ — 2U^(w2/2 -I- iuj, iv2,ioJ + ^2/2; ky, k, ki,by,b, hi, b') 

X v^{iuj + W2/2, iv2, -W2/2 - iuj] ky -ki+k, k2',k; b', b2',h, 62) 
-I- V^(iv 2 l 2 -I- iuo, iv2, iuj + 11^2/2; ky, k, ki;by,b, bi,b') 

X V^iioj + iv2l2,iuj + W 2 / 2 , -W 2 ] ky -ki+k, A;2', ^ 2 ; b', &2', &2, b) 
+ + W 2 / 2 , iw + W 2 / 2 , iv2] k, ky, ky b, by, bi, b') 

X + W2/2, iv2, —W2/2 — ioj-, ky — ki+k, k2yk; b', 62', ^2)] 

^CTph(»^3;fci',fc2',fci;&i',&2',&i,^2) ^ D^{iuj,k,b)D^{ivo+ii^,k + k 2 '-ki,b') 

b,b',k 


X V^{iuj + ivol2, iui + ivo/2, iv^-, k, k2',ki;b, &2', &i, &0 
X V^{iuj + ivol2, — iuj, ii/ 3 ; ky,k + ^ 2 ' — ki-,by,b', b, 62 ) 


( 6 ) 


( 7 ) 


( 8 ) 
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In this approximation, we can also write the full coupling 
function in the target bands at the end of the flow as 
follows, 


iv2, ^ 3 ; ki' , * 2 ', ki ; br , & 2 ', , & 2 ) 

— target—bands 

b'j ^2', b; &!', ^2', ^1, ^ 2 ) 

(*^2 ; kii , k 2 >, ki ; bii , &2' , , ^ 2 ) 

+‘J’CTph(*^3; kv,k2>,ki-,bii,b2',bi,b2) ■ ( 9 ) 

In the cRPA, the second and the fourth term on the right 
hand side are absent. Note that otherwise the flow for 
any A couples the different contributions $pp, ‘b^ph and 


B. Frequency structure of the vertex function 

As can be seen from the flow equations, the resulting 
vertex will depend on three frequencies. Unless one can 
insert this function directly in a many-body low-energy 
solver like CT-QMC, this rather rich frequency struc- 

I 


ture must be interpreted in some way. We propose to do 
this by separating the different channels and interpreting 
these as charge, spin and pairing terms of the effective 
interaction. In the following we explain how this decom¬ 
position is organized. This is done most easily in real 
space on the lattice of the effective low-energy model. 
As for the problems studied explicitly in the remainder 
of the paper, the effective model with contain only one 
orbital per site, we can henceforth ignore the band or 
orbital indices in the notation. 

Let us describe the fermions in the effective model with 
Grassmann fields cl ^{t) and Ci^o-('r) in imaginary time. 
Then the Heisenberg density operator at site i is given 
by 


p*b) = G>b)ci.<T(T), (10) 

a 

Pi{iv) = / dr exp(wr)pi(r) = ^ V' c\^{iui)cia{iui + iu). 

•^0 P 

( 11 ) 


A retarded density-density interaction term between site 
pairs (i,j) (symmetric in these indices, and each pair 
(i,j) counted once) can then be written as 




charge 


^ X] /charge(iGbj)P»(w)Pj (-*«") 


( 12 ) 


1 


1 

¥ 


-iv)cia{iuj + iv) ( 13 ) 


E EE4., ^{iujy)cl^^ y (iw2')14harge(*wi/, in; *W2', *2'|*wi, ii, za;2, i2)ci^,a' (ia;2)cii,cr(iwi) (14) 

cr,cr' 

iuJl,iuJ2 ^1^2 


with the general spin-rotationally invariant conpling function 

hcharge ; ^1' ? ' 1 ^ 2 ' ^^ 2 ' j ^^ 2 1^2) — ^ ^ /charge(^^7 ; '^2) ^ -^1 ^iu )2 —ii^,iuj 2 / ^13/ ,22 ■ 

iv 


( 15 ) 


Note that the conservation of Matsubara frequencies 
is implied by the Kronecker-ds in the last expression. 
Hence, for a given pair of sites (11,12), a density-density 
interaction with frequency transfer iv of the type (na) 
causes a fixed value /charge(iJ^, i,j) of the coupling func¬ 
tion V'charge(iwi',iiqia;2',i2'|*‘^i,U,iw2,12) for all exter¬ 
nal quantum numbers where the first frequency trans¬ 
fer Mandelstam variable, iuii — iwi', is the same. This 
transfer is the one between those incoming and outgo¬ 
ing frequencies of the legs that need to have the same 
spin index. Reversely, if we have a general effective cou¬ 
pling function V (fwi',Uq 1022',*1,*a;2, 12), e.g. ob¬ 
tained by integrating out some bands, we can average 
over the two other Mandelstam frequencies ioji — iuj2' 


and ioji -I- iuj2, or keep these fixed at specific values, in 
order to obtain an interaction that we can interpret as 
being of density-density type. Looking at the decom¬ 
position of the effective interaction in Eq. we could 
straightaway interpret (iiy; ky , ^2', b; &i', ^2', &i, ^2) 

as density-density contribution, but we will see that a 
part of this is also needed for the spin-spin part. Later, 
in a subsequent simplification, we will chose to ignore the 
fz/-dependence of that term as well, for being able to use 
it directly in an interaction Hamiltonian. 


We can proceed in a similar way with the spin part. In 
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imaginary time the Heisenberg spin operator is given by 


Si{iv) = / dr exp{ii'T)Si{T) 

Jo 


= ^ XI X! + w) 


( 17 ) 


Siir) 


1 

2 


E 




A retarded spin-spin interaction term between site- 
( 16 ) centered spins localized at i and j (again each pair 
counted once) can be written as 

I 


'S'spin — ^ 'y ( /spin(*^': b j)*S'i(lZ/) ■ Sj{ iv) 

= E /spin(w,bj)[^cJ^(iw)c]^,(iw')ci^/(ia; + w)cj>(ia;'-w) 


( 18 ) 


+ -cJg.(iw)cl^,(iw')cjCT'(ia;' - iv)cia{ioj + iv)] 


,12/,ti ,12 
<T,a' 


now with the coupling function 


14pin(*Wl',*l';*W2',l2'|jWi,Zi,iW2,*2) = “ E /spin 




( 19 ) 


= 4 E E ^ (ja;i/)c|^, (*a;2' )14pin(*wi/, H'; *^2', *2' |*wi, H, iu}2, 12)012 ,<t' (*W2)cii ,cr {iuJi) (20) 


“f ^ ^iuj2 —112,1(^2* 


( 21 ) 


Therefore, the spin part involves a dependence on both 
Mandelstam frequency transfers. For fixed iv we can 
obtain a relation between the terms in the decomposition 
® and /charge/spin(iJ^,b j) as follows, 

^dfTh (*i^;HA2,'il,i2) = /charge(w,H,*2) “ ^spin(w, *1, 12 ) 

( 22 ) 

^crph(*i^;*2,n,'il,i2) = -x/spin(w,Ji,i2) 

( 23 ) 

This allows us to understand the data from the 
fRG from integrating out the high energy bands, 
given by the functions *2, U, *2) and 

‘l’crph(*^> *2)*1)*1) *2) as density-density- and spin-spin 
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interacitons /charge/spinb i)- These functions have 
to be added to the bare direct interaction Ttarget-bands 
between the low-energy degrees of freedom. 

The interpretation becomes simplified if one throws out 
the frequency dependence of the effective interaction 
functions /charge/spin (*J^, b j)- Then one can formulate an 
effective interaction Hamiltonian. For comparison, let us 
consider an interaction Hamiltonian of the form 

H ^'^UijPiPj + Sj ( 24 ) 

(id) (id) 

with instantaneous density-density and spin-spin interac¬ 
tions. The corresponding coupling function is, after in¬ 
serting Grassmann fields for the equal-time density-and 
spin-operators, given by 


F(icdi/, H' 1 ^( 3 ^ 2 b ^2' 5 H 5 ^^2, ^2) — y ( ^Gq, ^'I2,/jb/ diuj 2 —ii 2 ,iLJ 2 i di'^i 2 d 2 f 

L \ / 


2 'li,/ ^i^l+ii2d^^2' ^i^2 —ivd^V ^ild2' dv 


( 25 ) 


We can see that this has the same structure as the sum 


of Eqs. [T 5 ]and[ 2 T] One way to proceed is now to identify 

Uij = /charge(w = 0, i,j) , and ( 26 ) 

Jij — /spin(^^' — 0 ,bj) ■ ( 27 ) 
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This will most likely overestimate the effects of the small 
^-contributions, but will still allow to obtain qualitative 
answers about the tendencies caused by the renormal¬ 
izations that affect /charge/spinbj)- Alternatively one 
could use w-averages of the functions /charge/spin(*^^1 hj)- 


We note that so far we have restricted the study to local 
charge and bilinear spin. It is however straightforward to 
generalize the considerations and identifications above to 
nonlinear charge and spin fields that live, e.g., on sites 


Then we obtain the relations 

^dtTh *lb *2', n, * 2 ) = /charge(w, *1', *2', *1, * 2 ) 

^y'spin f 2q f 15 f 2) 

( 29 ) 

The identifications are still taken with respect to the fre¬ 
quency structures. 

Finally we can also consider pair-pair interaction of the 
type 


‘S'pairing — 03 ^ ^ ^ ^ /pairing(f fl') f 2') *1) f 2 )Cj^, g. (iW ^/( iu! il')Ci2a'{ iUJ ( 30 ) 

,i 2 / 41^2 

a ,<j' 

^ ^ X] X] X~! .ct' (»^2')f^pairing(f^l', fid »^2', ^2' I 1 1 , f (^2, (31) 

(J (yl 

*1.^2 


with the general spin-rotationally invariant pairing interaction 

hpairing('^^l^ ? ; ^^2'; ^‘ 1 ’ 5 5 ^ 2 ) — E /pairing (iz2, fi', ^2') *1, ^ 2 ) <5 (32) 

iv 


In the employed vertex decomposition we get 
^pir^(*'^;*lb*2 ',*l,*2 ) = /pairing(w, *1 ',*2',*1,*2 ) • ( 33 ) 

This completes the representation of the effective inter¬ 
action in terms of charge-charge, spin-spin and pairing 
interactions. 

III. RESULTS ON MODEL SYSTEMS 

Here we describe the numerical results obtained by ap¬ 
plying the cfRG scheme described so far to two different 
few-band models. In the first model, we find only minor 
deviations from the cRPA results, while in the second 
example, strong changes occur. 

A. Graphene model 

1. Hamiltonian 

As a first application we discuss in the following a multi¬ 
band model of graphene<^ The unit-cell of a graphene 
layer is shown in Fig. [S] It consists of two sites with each 
3 sp^-hybridized orbitals, forming the cr-bonds. In addi¬ 
tion there is one Pz-orbital per lattice site that undergoes 
TT-bonding with the Pz-orbitals on neighboring sites. The 


tj-bands shall serve as high-energy bands to be integrated 
out, while the 7r-bands are left for the LEEM. These two 
target bands also feature the well known massless Dirac 
spectrum at low energies. 



FIG. 2. (color online) sp^ hybridized orbitals in graphene on 
neighboring C sites. We also indicate the hopping parameters 
included in the Hamiltonian. 

The Hamiltonian of our model consists of a free part Hq 
and an interaction part Hint 

= iJo + Hint ( 34 ) 

The free part is given by 








X! 4,a,0^fe,o,0 + X! X! X! X! (4,a,z'^fc,oj + ^■'^•) 


3 =A,B 


+ (4A,0^fe,B.O + ^-C-) + 


a=A,B i^O 
ik&i t 


a=A,B O^i^j^O 


e ^C 


JkS2J 


k,A,l'"k,B,l ^ ^^k,A,2''k,B,2^''k,A,3''k,B,3 


+ . „cn 


h.c. 


with 

5i = (l,0), ^2 = ^(-l,\/3), (?3 = i(-l,-\/3) 

7 £= E 

i^l, 2,3 

l 7 fcl = Y 3 + 2 cos - (52)^ + 2 cos - 83)^ + 2 cos [k{62 - 63)^ 


( 35 ) 

( 36 ) 

( 37 ) 

( 38 ) 


Here, we choose to compute the values of the parameters 
above in a Slater-Koster-formalism (see, e.g. Ref. IH), 
with localized hydrogen wave-function parametrized by 
an adjustable Bohr radius, as in RefUf]. The size of the 
matrix-elements is given by^i 


= -ll.OTeR ( 39 ) 

e,, = - 13 . 84 ey ( 40 ) 

iintra = - 2 . 776 ^ ( 41 ) 

= - 2 AeV ( 42 ) 

= -l 2 MeV ( 43 ) 


The resulting band structure with two 7r-bands at the 
Fermi-energy and six high energy cr-bands is shown in 
Fig. [ 3 ] The F, M and K-points in the first Brillouin-zone 
are shown in Figure [H 



FIG. 4. The first Brillouin zone of hexagonal lattice for the 
graphene model. 



FIG. 3. (color online) The band-structure of graphene corre¬ 
sponding to the noninteracting Hamiltonian (I35F There are 
two TT-bands at the Fermi-energy (red) and six high energy 
CT-bands (blue). 


Note that the matrix elements between tt- and cr-orbitals 
are equal to zero. This is due to the mirror symmetry 
Mz or z-reflection symmetry that takes the coordinate 
perpendicular to the graphene layer, z, to —z. The <j or¬ 
bitals are even under this symmetry while the rr-orbitals 
are odd. By a nonzero perpendicular electric field E or 
by deposition of the graphene sheet on a substrate, the 
Mz symmetry can be broken explicitly. This leads to an 
additional term in the Hamiltonian of the form 


Httct — Ittct E E (4,a,i'^fe.a.O + ^.c}j ( 44 ) 

^ a,i^Q 
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The interaction part is taken as 


i?int=E E 


n a—A,B 


tEr7r^n,a, 0^n,a,0 “I” Uaa E ^n,a,i^n,a,2 


H“ ^CTTr ^ ^ ^n,a,oOn,a,z 

“t“ t^intracr (‘^n,a,l‘^n,a,2 H“ ^n,a,l^n,a,3 “t” ^n,o,2^n,a,3) 

(45) 


The interaction matrix elements are given by {V (r —r^) = 
Bohr-radius oq = 5.29 10“^^m). 


b^TTTT — 

Uaiz — 

^intracr — 


J d^r J (i^r^^/!o(r)^r(r — F')^o(^)^ ~ O.I 97 E 
J d^r J d^r''il}i{f)^V{f— r')'0i(r')^ Ri 0.204^ (j ^ q) 

J J d^f^'ipi{f^‘^V{f — r')il^o{r')‘^ 0.171— (i ^ 0) 

f d^r f dV?/'i/2/3(^)^r(f- i^)V'2/3/i(^)^ ~ o.iseE 


(46) 

(47) 

(48) 

(49) 


with ^ = 27.2ey and Z* = 3.25. Ut^tt is in good agree¬ 
ment with Ref. m. 

In principle, the general representation of the Coulomb 
interaction in the localized-state basis gives rise to many 
other terms, e.g. such terms with wave function centered 
at more than two sites. We ignore all other terms that 
are not of density-density- and two-center-type, as the 
values for these integrals turn out to be small. Note 
however that this statement hold only on the ’tree level’ 
of the bare integrals without quantum fluctuations. One 
of the main points of this paper is that loop-corrections to 
these parameters can actually generate noticeable terms 
that are not of density-density type. 


2. Effective interaetions in the ir-seetor 


(a.) 



77 77 


(b.) 



CT CT 


(c.) 



CT a 




Urr 


Uaa ) f^ntracr 


FIG. 5. (color online) The bare one- or two-center interaction 
vertices taken as nonzero. The spin projection is conserved 
along the short side of the rectangle. The interactions can be 
classified in three types: a) Terms Ut^tt which couple different 
TT-orbitals, b) the terms oc Uan which couple a- and tt- orbitals 
and c) the terms oc Uaa,Uintracr which couple only a-orbitals. 


Now our aim is to derive effective interactions between 
the TT- orbitals by integrating out the high energy cr- 
orbitals. This integration can be done either by using 
the constrained RPA method or the wick-ordered fRG. 
The question of our concern is if the latter leads to a 
qualitatively new picture, i.e. if non-RPA terms are im¬ 
portant. 

As argued above the tt- and the u-orbitals are not cou¬ 
pled in the free part Hq as long as the mirror symmetry 
Mz holds. This symmetry decoupling has profound con¬ 
sequences on the perturbation series that is summed in 
integrating out the cr-bands. In Fig. [S]we show the bare 
interaction terms (I46M9I) , which can be classified in three 
types: (a.) the terms oc which couple 7r-orbitals, (b.) 
the terms oc Ua-K which couple a- and tt- orbitals and (c.) 
the terms oc Uaa-, lAntracr which couple only cr-orbitals. 


Now we consider the effective interactions in the 7r-sector 
in second order perturbation theory, with the internal 
lines constrained to contain at least one propagator in 
CT-sector. In Fig. |6] a) and b) we show possible contribu¬ 
tions to the effective interaction from the particle-particle 
diagram. As the u-orbitals are integrated out at least one 
inner line of this diagram has to belong to the tr-sector. 
But as shown in Fig. [Sj the bare interaction vertices 
corresponding to the Hamiltonian (HSl) have a different 
structure than the allowed interaction vertices in Fig. [51 
Therefore we do not get any nonzero contributions to the 
effective TTTr-interaction from the particle-particle channel 
in second order perturbation theory. A completely anal¬ 
ogous argumentation holds for the crossed-particle hole 
channel and the last two ’vertex correction’ diagrams of 
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the direct particle-hole channel (cf. Fig. [T]). Of course 
the vanishing of this contributions depends on neglecting 
bare interactions other than of density-density type, but 
if these were nonzero, only minor corrections would oc¬ 
cur. In principle, processes within the high-energy bands 
can also these non-density-density terms, as will also be 
visible in our data. But these terms just lead to small cor¬ 
rections, as their magnitude is controlled by the energy 
gap AE between the target bands and the high-energy 
bands. 



TT a 7T 7T TT 7T 


FORBIDDEN FORBIDDEN 



FIG. 6. (color online) Different diagrams of second order as 
appearing on the right hand side of the fRG equation for the 
effective interactions. Due to the symmetry-restricted set of 
initial interactions, the only nonzero contribution with 4 ex¬ 
ternal TT-band legs comes form the RPA-type diagram in (c). 
The contributions from the particle-particle diagram, from 
the crossed particle hole diagram and from the vertex correc¬ 
tion diagrams in the direct particle hole channel cancel. 

Thus it turns out that the only significant renormaliza¬ 
tion that is generated in second order perturbation theory 
with the bare interaction vertices from Hamiltonian (I45|) 
results from the first direct particle-hole diagram shown 
in Fig. ^c). This diagram is of course also included in 
the cRPA approach, i.e. on this level of approximation 
both approaches are equivalent. 

Beyond second order perturbation theory, the picture be¬ 
comes slightly more difficult. For example interactions 
exclusively in the cr-sector that can be generated from all 
channels can be inserted into the inner part of the dia¬ 
gram of Fig. ini c. However, the numerical analysis for the 
realistic parameters listed above does not reveal sizable 
deviations from the cRPA. 

In Fig. [7] we show data for the effective interactions in 
the target bands, for the case where all external legs are 
chosen at one AT-point in the same 7r-band. It can be 
clearly seen that the differences of the cfRG to the cRPA 
are at best of minor quantitative type. In the cRPA cal¬ 
culation the vertex only depends on the second transfer 
frequency iv 2 corresponding to the frequency dependence 


generated by the first direct particle-hole bubble. In the 
cfRG calculation no additional frequency dependence on 
ivi or is generated and a peaked structure depend¬ 
ing on iv 2 is observed. The structure is qualitatively the 
same as in the cRPA calculation. 

Of course, this is only a statement about a single point 
in momentum space, but it illustrates the diagrammatic 
argument that non-cRPA corrections are weak. Thus we 
can conclude that also in more ab-initio based studies 
for this system like done in Ref. [s^ the result of the 
cRPA-analysis should come very close to the answer one 
would obtain when more than the cRPA diagrams are 
considered. 

This picture can be expected to change parametrically 
when we include a symmetry-breaking term ()44p . Then, 
the high-energy bands will obtain admixtures of the tt- 
orbitals and vice versa. This opens new interaction chan¬ 
nels, as now the formerly forbidden terms lead to contri¬ 
butions parametrically large or small in the strength of 
the symmetry-breaking. This situation will occur e.g. in 
layered graphene, e.g. in bilayers, or for the outer lay¬ 
ers or asymmetric bands of trilayer systems. Similarly, 
putting monolayers on a substrate or electric field in z- 
direction will violate the z-reflection symmetry as well. 
We here try to capture these effects qualitatively by in¬ 
troducing a nonzero hopping that now couples the 
two sets of bands. 

In Fig. |S] we compare the fRG and the cRPA cal¬ 
culation of the coupling function component Vmi = 
V[K,tta\K,tta\K,'ka]K,'ka) in the target 7r-bands at 
zero frequency transfer as function of t,rCT- K is one of the 
Dirac-Points shown in Fig. |4l The bare matrix element 
= 17.4eP becomes screened both in the cRPA and 
in the fRG calculation. For t,r(T = 0 the screening effect 
is larger in the fRG calculation than in cRPA, which is 
caused by the corrections beyond second order pertur¬ 
bation theory. For 0 a different behavior in both 

approximations is observed. In the cRPA calculation, the 
screening becomes larger when is increased, while in 
the fRG calculation we obtain an anti-screening of the in¬ 
teractions which is mainly caused by the crossed particle- 
hole channel (cf. the channel-decomposed screening be¬ 
havior in Fig. [9] that shows the increase in the crossed 
particle-hole channel coincident with less reduction in the 
direct particle-hole channel). Note however that this dif¬ 
ferent behavior with becomes significant at values 
of t,rCT that are very large compared to the bare energy 
scales of the Hamiltonian. Furthermore, at t,r(T ^ 0, the 
corrections to cRPA make the cfRG-result more similar 
to the cRPA data. So, in the present example, the cRPA 
results appear to be quite robust even if the two sets of 
bands are allowed to couple. 

An important learning from this study is however that 
the symmetry of a problem can decide if there are cor¬ 
rections to cRPA or not. If the z-reflection symmetry is 
broken, new terms enter the perturbation series for the 
effective interaction which contain physics that go be¬ 
yond screening arguments. In principle this opens a new 
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FIG. 7. (color online) The vertex V{K,nA', KjTyaIKjTya’, KjTva) — U in cfRG (left) and cRPA (right). 




FIG. 8. (color online) Gomparison of the fRG and 
the cRPA calculation of the matrix-element Vim = 
V{K ,'KA', K ,nA\K,nA', K ,'ka) at zero frequency transfer as 
function of the hopping 


FIG. 9. (color online) Comparison of the three interac¬ 
tion channels = V{K ,tva', K,tva\K,tva’, K ,'jta) at zero fre¬ 
quency transfer for i = pp, dph, crph with the cRPA calcula¬ 
tion. 


variant of tuning the interaction parameters of a two- 
dimensional system, although the precise impact of the 
symmetry breaking is possibly not intuitively clear. 


B. One-dimensional chain 

1. Hamiltonian 

As a second example we consider a simplistic model for 
one-dimensional chain, inspired by copper oxides. It con¬ 
sists of ’copper’-centered 4s- and 3d-orbitals and ’oxygen’ 
sites with 2p-orbitals extending towards the copper sites 
(cf. Fig. [TUll . This model can be considered as more 
generic than the graphene example as there is no symme¬ 


try difference between the bands near and further away 
form the Fermi level. 

The Hamiltonian is given by 

H = Ho + ffint (50) 

Ho = esJ2 Pk,crPk,a 

k,<7 k,(7 k,<7 

“h tsp E [1 -I- exp(ffc)] (^sl^^pk,cr + H.C.^ 

fc,(T 

^dp E [1 -I- exp(iA:)] (dl^^pk,a + H.c^ (51) 

k,(T 
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Hirit — Us + Ud dl.^di^-fdl^^di^i 

i i 

+ Up^pl^p^,^pl^p,^l + Usd ^ sl^S^,ad\,^,di^a' 

i i,(T,cr' 

(52) 

We choose e* = 2, = -0.5, Cp = -2, tsp = 1, tdp = 1, 

Us = Up = Ud = U = 1.5, Usd = U/3, all measured 
the same energy unit. The temperature is chosen as /3 = 
1/T = 5. 

The hopping element between the copper 4s- and 3d- 
orbital is equal to zero, due to symmetry. The spectrum 
of the free part Hq is shown in Fig. [TTJ It consists of two 
high-energy bands above and below the Fermi level and 
one conduction band. They are called s-like, d-like and 
p-like according to their main orbital contribution. With 
parameters listed above, the d-like is (roughly) half filled. 



Cu O Cu O 

FIG. 10. (color online) The one-dimensional chain consists of 
’copper’ sites with 4s- and 3d-orbitals and ’oxygen’ sites with 
2p-orbitals. The hopping matrix element between the s- and 
d-orbitals on the same site is equal to zero due to symmetry. 

In the following we integrate out the p-like and the s- 
like band to get an effective theory for the d-like band. 
Note that now, due to the overlaps tdp between the d- 
orbitals and the p-orbitals, the d-orbital also has some 
weight in the high-energy bands. This opens the way for 
non-RPA contributions in the integration over the high- 
energy bands. Hence we expect stronger deviation of the 
cfRG from the cRPA. 


8. Effective Interactions 

As the model is one-dimensional, we can also take first 
steps to explore the wavevector dependence of the effec¬ 
tive interactions. We do this by dividing up the one¬ 
dimensional Brillouin zone in N patches. Here we use 
N = 2 and N = 4, but with slightly higher numerical 
ambitions, this number can be increased. To be more 
precise, in a Scheme 1 we use two fc-points at fc = 0 and 
k = n and 41 Matsubara frequencies. In a Scheme 2 we 
use the fc-points k = —7r/2, k = 0, k = 7r/2 and k = tt 
and 11 Matsubara frequencies. We restricted the internal 



FIG. 11. (color online) The spectrum of the free part Hq of 
the Hamiltonian (|50l). It consists of two valence bands and 
one conduction band. They are called s-like, d-like and p-like 
according to their main orbital contributions. 


fc-sums in the flow equations also to this discrete set of 
fc-points, which is a valid approximation due to the rela¬ 
tively weak dispersion of the bands compared to the band 
gaps. So, technically, our calculation corresponds to a fi¬ 
nite size cluster and by Fourier transformation we can 
obtain local and nonlocal effects up to three neighbored 
lattice sites. 

In Fig. [12] we show the three interaction channels 
i,j, k,l]d, d, d, d) with i = pp, dph, crph in posi¬ 
tion space, comparing cfRG in the left plots with cRPA 
in the right plots. Of course, in the cRPA only the direct 
particle-hole channel contributes by construction. Note 
that $crph(w = 0) is positive in all cases in the cfRG, 
which corresponds with Ea. (l?51) to a negative (i.e. fer¬ 
romagnetic) exchange interaction J. Below we will show 
that this generated exchange interaction leads to a polar¬ 
ized ground state at least in a small system. In the cRPA 
calculation one has 4>crph(*i^) = 0 and therefore the ef¬ 
fective J according to our decomposition of the effective 
interaction is zero. The nearest-neighbor direct particle- 
hole term in the second row has a distinct behavior as a 
function of the transfer frequency due to a sign change 
between the fc = 0 and the fc = tt term. This is already 
visible in second order perturbation theory. 

In Fig. [T3| we show contour plots of the frequency- 
dependent part of the full vertex as sum of the three 
channels in the W 2 -ir' 3 -plane (with ivi = 0) for the cRPA 
and the cfRG calculation. Here the differences can be 
clearly read off. In the cRPA calculation the vertex only 
depends on transfer frequency W 2 , while in the cfRG cal¬ 
culation more frequency structure evolves, which leads 
to additional peaks at ^2 (and also ivi, not visible in 
the Figure). These additional features are often compa¬ 
rable or larger than the structures that are found in the 
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FIG. 12. (color online) Effective d-band interactions at [/ = 3, in units of tsp- First row: The three interaction channels 
1, 1,1,1; d, d, d, d) with i = pp, dph, crph in £RG (left) and cRPA (right). The indices T, 1,1,1’ refer to lattice sites, so 
this is the local interaction. Second row: The three interaction channels {iv, 1, 2,1, 2; d, d, d, d) with i = pp, dph, crph in fRG 
(left) and cRPA (right). The site combination T, 2,1, 2’ is such that the spin projections on site 1 and 2 remain both unchanged 
in the interaction process. Third row: The three interaction channels 1, 2, 2,1; d, d, d, d) with i = pp, dph, crph in fRG 

(left) and cRPA (right). The indices 1 and 2 refer to neighboring lattice sites. The site combination T, 2,1, 2’ is such that the 
spin projections on site 1 and 2 are exchanged in the interaction process if the incoming spin projection is zero. 


cRPA calculation. In particular, the cfRG effective inter- C. Impact of new terms in effective Hamiltonian 
actions can also enhance the bar interactions and cause and interpretation 

’anti screening’, while the cRPA screening always tends 
to suppress the interactions. 

We can now ask what the consequences of the additional 
interactions in the LEEM beyond the density-density 
terms are. One way to proceed with the full frequency 
structure would be to use CT-QMC to solve a small clus- 
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FIG. 13. (color online) Effective interactions at 1/ = 3. First row: The vertex = Q,iu 2 ,wz', 1,1,1, l;d, d, d, d) — [/ in 

clRG (left) and cRPA (right). Second row: The vertex = 0, ii/2,ii'3; 1, 2,1, 2; d, d, d,d) — U in cfRG (left) and cRPA 

(right) Third row: The vertex = 0, iv 2 ,iv 3 \ 1, 2, 2,1; d, d, d,d) — U in cfRG (left) and cRPA (right). 


ter (that could also be embedded via DMFT). As we 
do not have this heavy machinery available here, we re¬ 
sort to solving an approximate Hamiltonian that contains 
the new terms in exact diagonalization, again on a small 
number of N sites corresponding to the wave vector dis¬ 


cretization described before. 

In order to solve the effective theory for the d-like band 
we construct an effective Hamiltonian with the vertex 
functions at zero frequency transfer as interaction matrix 
elements, along the lines of Sec. Ill HI More precisely we 
write 


H = 






E 




(53) 
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where the interaction constants are given by 


^{ivi = 0,iv2 = 0,ii/3 = Q.iv,i2’.,11,12, d,d,d,d) 


(54) 


These interaction terms can be interpretated in terms of 
charge, spin and pairing interactions according to Eqs. 
I28I29I and [33l Instead of using the interaction terms at 
ivj = 0, we could also use Wj-averages in each direction 
j = pp, dph, crph which would lead to slightly smaller 
interaction terms. But on a qualitative level we do not 
expect differences in the results and the detailed choice of 
approximation will be of minor importance. In this way 
all retardation effects are neglected, but we are able to 
solve the effective Hamiltonian exactly. In both schemes, 
i.e. for TV = 2 and for iV = 4, we found that for a high 
enough interaction strength U the groundstate at half 
filling of the remaining d-band is given by the respective 
highest spin state. In scheme I with 2 sites the ground- 
state spin is given by S' = I and in scheme 2 with 4 sites 
it is S = 2. As a measure for the strength of the spin- 
spin interaction we used the ’spin gap’ A between the 
groundstate and the first excited state (with spin S = 0 
in scheme I and S = I in scheme 2). In Fig.[T4]we show A 
as function of U in both schemes for three different types 
of approximation. In the simplest ’tree-level’ approxi¬ 
mation we used only the bare interaction (i.e. without 
integrating out the excited energy levels) in the effective 
Hamiltonian for the exact diagonalization. In the more 
advanced approximations we integrated out the excited 
levels by the cRPA and the fRG. 

It turns out that the high-spin groundstate is already 
found in the tree-level approximation, where we use the 
bare interaction parameters for the effective model. So 
the tendency toward a high-spin state is not purely an 
effect from integrating out the high energy s- and p- 
like bands. Note however that this part of the interac¬ 
tion is frequency independent and cannot be interpreted 
straightforwardly as an interaction term like Eg. (11811 . 

Compared to this bare case the energy gap decreases in 
both discretization schemes (TV = 2 and TV = 4) when we 
integrate out the high energy levels in the cRPA scheme 
and increases in the cfRG scheme. In all cases there is 
a minimal threshold value of U above which A becomes 
zero (and negative for smaller values of t/). In scheme 1 
with just two wave vector states this can be easily under¬ 
stood, as there are only two energy levels in the effective 
Hamiltonian at fc = 0 and fc = tt with different onsite- 
energies e/c= 7 r < efe=o- For a small interaction strength 
both electrons occupy the k = 7r-level forming a spin 
singlet state. For a larger interaction strength the dou¬ 
ble occupancy becomes energetically unfavorable and the 
electrons start to occupy the k = 0-level. Due to a posi¬ 
tive exchange interaction between both levels they form 
a high spin state. 

Very roughly, this increased tendency toward spin 
polarization can be understood as an inter band 




FIG. 14. (color online) Energy gap A between the high-spin 
groundstate and the first excited state as function of U in 
scheme 1 (above) and scheme 2 (below) using three different 
approximations for the interaction terms in the exact solution 
of the effective d-band-model. In the simplest case we used 
just the bare interaction as input for the effective Hamilto¬ 
nian. This is compared to the energy gap from the effective 
Hamiltonian with interaction parameters from the cRPA and 
the fRG calculation. 


RKKY(Ruderman-Kittel-Kasuya-Yosida )-mechanismi^ 
In our case the role of the conduction band that medi¬ 
ates the spin-spin interaction is taken by the high-energy 
bands, while the localized moments that interact are rep¬ 
resent by the (still mobile) electrons in the LEEM target 
band. In the usual RKKY effect for a metal, the spin 
susceptibility of the metal causes a power-law spin-spin 
interaction with a sign structure controlled by the Fermi 
surfaces geometry. Only for distances less than the Fermi 
wavelength, the mediated interaction is ferromagnetic. 
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beyond this scale the sign of the interaction oscillates. 
In our case, there are no particle-hole pairs that happen 
solely near the Fermi level, and the wave vector structure 
of the particle-hole bubbles is quite mild. Hence there are 
no sign oscillations, and at least within the range of our 
calculation, the interband-mediated spin-spin interaction 
remains ferromagnetic. 

Note that we have here studied a weakly-correlated elec¬ 
tron system and ignored some known effects. For exam¬ 
ple, a one-dimensional half-filled conduction band with 
repulsive interactions will become insulating due to Umk- 
lapp processes. This effect is not quite visible in our solu¬ 
tion of the small cluster which is dominated by finite-size 
gaps. For larger (for Cu more realistic) local interac¬ 
tion parameters the Mott-insulating state at half filling 
will feature antiferromagnetic exchange interactions be¬ 
tween the localized electrons. This mechanism will very 
likely outweigh the ferromagnetic tendencies described 
above. Yet, on a quantitative level, the corrections to 
cRPA should remain present. 


IV. CONCLUSIONS 

We have studied the perturbative calculation of effective 
interaction parameters in multiband models with short- 
ranged interactions for electrons in solids. We have de¬ 
scribed how the functional renormalization group (fRG) 
can be employed to obtain these interaction terms of 
the low-energy effective model (TEEM) and how this 
treatment goes beyond the commonly used constrained 
random phase approximation (cRPA). Eor two models, 
we compared the cfRG scheme and the cRPA. In the 
graphene-type model, both approximations gave similar 
results. The reason for the absence of noticeable correc¬ 
tions to the cRPA could be identified as the different sym¬ 
metry with respect to z-inversion between the 7r-orbitals 
forming the conduction band and the cr-orbitals forming 
the bands that are integrated out. When the z-inversion 
is destroyed by additional terms in the free Hamiltonian, 
new terms beyond cRPA enter the cfRG flow and the 
LEEM interactions from the cfRG differed considerably 
form the cRPA. This opens a new way how effective inter¬ 
action parameters of a system can be tuned, i.e. when the 
embedding of the system breaks a certain symmetry and 
hence makes non-cRPA contributions nonzero. Rather 
generally, there are two other important channels beyond 
the cRPA in this setup: the particle-particle channel that 
in general suppresses the interactions further, and also 
the crossed particle-hole channel that is responsible for 
magnetic tendencies and that works to increase the effec¬ 
tive interactions. Depending on the model parameters, 
not much of the original cRPA screening may be left in 
the cfRG effective interactions. 

There are various paths that need to be pursued now: 

1) The solution of the LEEM should be improved, al¬ 
ready for the rather simplistic full models analyzed here. 
Then the meaning and impact of the terms beyond cRPA 


will become more clear. Besides employing state-of-the- 
art impurity solvers for clusters together with embedding 
schemes, nanosystems may represent a fruitful test arena, 
as additional benchmarking should be possible. 

2) The formalism presented here needs to be cast in the 
more general framework of ab-initio calculations. Most 
likely, it will not be possible to capture all bands in the 
cfRG, as this would lead to an extreme numerical load, 
but ist is well conceivable that improved approximation 
scheme will treat the high-engery bands closest to the 
target bands in cfRG, while for higher-lying bands cRPA 
or simpler approximations are possible. Note also, that 
the inversion of the dielectric constant matrix that makes 
the cRPA potentially a numerically complicated task is 
not necessary in the cfRG. Here the solution of the dif¬ 
ferential equation takes over this part. 

3) We have focussed on the effective interaction of the 
low-energy model. Of course, self-energy effects might 
also be an important factor. In principle, the flow of the 
self-energy is also described by the fRG formalism, we 
have however not yet analyzed how this should be imple¬ 
mented in the question of deriving effective low-energy 
models. A recent study of Shinaoka et ali^ studies ex¬ 
actly this issue by powerful (extended) DMFT methods. 
The main upshot of their work is that a self consistent 
treatment of the energetics of the bands is essential, and 
that interacting propagators should be used if one tries 
to compute the effective interactions. It will be interest¬ 
ing to see how our work and the approach by Shinaoka 
et al. can be tied into a comprehensive picture. 

4) In the present work, the initial interactions were 
taken to be short-ranged while in the ab-initio context 
Goulomb interactions should be considered. One of the 
reasons for this simplification is that any additional 
wave vector dependence, in particular one related to 
longer-ranged interactions, increases the numerical 
effort. In this paper, our main goal was to expose 
that the inclusion of the additional diagrams can cause 
substantial changes, and this is clearly visible already 
in the frequency structure we focussed on. We however 
expect that also the wave vector structure of the effective 
interactions will turn out much richer in cfRG than in 
cRPA. This can already be inferred from computing 
second-order corrections. Another thing than can be 
expected is that in the limit of small-wavector transfers 
between legs 1 and 1', the deviations from cRPA will 
become smaller, as then the RPA diagrams of a given 
order will be helped by at least one factor 1 /qa (with 
lattice constant a) compared to any other diagrams. 
Note however that one of the main applications for the 
effective interactions is to distill from them effective local 
and short-ranged interactions. Hence, one will typically 
average over the wave vector transfers and specialties 
at small transfers will not drastically influence the 
results. Therefore the deviations from cRPA for generic 
momentum transfers will not be negligible. 
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